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We summarize the current status of our numerical results for the gluon and ghost propagators 
and for the Kugo-Ojima confinement parameter in quenched SU (3) lattice Landau gauge theory. 
The data for the propagators are compared to our results obtained in the case of full QCD, simu- 
lated using two flavors of dynamical clover-improved Wilson fermions. We demonstrate that the 
infrared behavior of the ghost propagator is consistent with the Kugo-Ojima confinement crite- 
rion. Explicit violation of reflection positivity by the gluon propagator is shown. Additionally, 
we present results of a running coupling constant both at low and at large momenta. 
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Lattice study of QCD Green's functions 



1. Introduction 

In the last years progress has been made to improve our understanding of QCD Green's func- 
tions in different gauges. This is important, because it might allow us someday to describe all 
hadronic features in the continuum directly in terms of those functions. A deeper knowledge of 
QCD Green's functions would help us also to understand the basics of nonperturbative phenomena 
like gluon and quark confinement or dynamical chiral symmetry breaking from first principles. For 
example, the realization of the Kugo-Ojima confinement scenario [[Tj] in QCD in covariant gauges 
is encoded in the infrared behavior of the gluon and ghost 2-point functions. Therefore, the in- 
vestigation of QCD Green's functions is not only of interest for a coherent description of hadronic 
states but also for an understanding of confinement. 

At large (Euclidean) momenta QCD Green's functions can be described in terms of perturba- 
tion theory. However, at intermediate and low momenta nonperturbative methods are indispensable 
to arrive at a complete picture. There are different approaches to deal with nonperturbative QCD, 
but they all have their own limitations. Therefore, a comparison of results is important. For exam- 
ple, studies of truncated Dyson-Schwinger equations (DSEs) for the gluon and ghost propagators in 
Landau gauge came to the conclusion that at very low momenta these propagators are governed by 
power laws with interrelated exponents which then results in a coupling constant running to a non- 
trivial infrared fixed point [§]. This infrared behavior has been confirmed independently by studies 
of exact renormalization group equations [Q] and investigations based on the Fokker-Planck-type 
diffusion equation of stochastic quantization [§]. At intermediate momenta, however, predictions 
e.g. based on DSE studies are uncertain due to the truncations involved, but in this region lattice 
Monte Carlo (MC) simulations can help. 

2. Lattice results for the gluon and ghost propagators 

We have performed MC simulations in the quenched approximation of SU (3) lattice gauge 
theory at four values of the coupling constant (j8 = 5.7, 5.8, 6.0, 6.2) using the standard gluonic 
Wilson action. Thereby, the lattice size ranged between 16 4 and 56 4 . To study also the influence 
of fermions we have analyzed gauge field configurations generated with Nf = 2 dynamical flavors 
of clover-improved Wilson fermions using the same gauge action 1 . All gauge configurations have 
been fixed to Landau gauge by maximizing the Landau gauge functional 

1 4 

F u[8} = T^LL^ Tl8U -^ with g U x ,»=g x U x , >1 gl +(l (2.1) 

x n=\ 

using either over-relaxation or Fourier-accelerated gauge-fixing. Of course, the maxima are not 
unique, but we will neglect the influence of Gribov copies in the following. The reader is referred 
to Ref. [|6j] for a more systematic account of the Gribov ambiguity in this context. 

After gauge-fixing we have calculated the gluon and ghost propagators on all (gauge-fixed) 
configurations. The definition of both Green's functions and their dressing functions is standard 
and our notation can be found in detail in Ref. n. To reduce finite volume and discretization 
effects we have applied cone and cylinder cuts to our data [0]. 

'Those configurations were provided to us by the QCDSF collaboration. 
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Considering first the quenched case, in Fig. [I] 
we show data for the ghost dressing function J(q 2 ] 
(renormalized at /I = 4 GeV) as a function of the 
momentum q 2 . As expected the ghost dressing 
function seems to diverge at vanishing momen- 
tum. We have tried to fit the infrared power law 
as expected from DSE studies to our data. How- 
ever, within the region of lower momenta a power 
law does not describe the data that well. The 
ghost dressing function seems to increase loga- 
rithmically as the momenta become small. In any 
case, the infrared exponent extracted from the fit, 
K « 0.2, is much smaller then it is expected to 
be, namely K 0.59 ||5|, ||]. With respect to our 
data for the gluon propagator, shown in Fig. ||, 
we think that we are still not in a region of mo- 
menta where the mentioned infrared power law can be verified. For this the gluon propagator D 
has to vanish at zero momentum which cannot be concluded from our current data (see Fig. |2|). To 
what extent finite volume effects can be blamed for such a (different) behavior found on the lattice 
needs to be clarified in future yet. 

To improve our understanding of lattice Landau gauge theory we have also studied the influ- 
ence of (clover-improved) Wilson fermions on the ghost and gluon propagators. We find that the 
ghost propagator stays almost unchanged if fermions are added to the gauge action. In contrast to 
this, fermions affect the gluon propagator at large and intermediate momenta, in particular where 
the gluon propagator exposes its characteristic enhancement compared to the free propagator. To 
illustrate this, in Fig. [3] we show our data for the gluon dressing function obtained for quenched and 
full QCD. 

As a part of the project we also performed simulations at j3 = 6.0 using large asymmetric 
lattice sizes, namely 16 3 x 128 and 24 3 x 128. By comparing the data obtained on these lattices for 



Figure 1: The ghost dressing function J for the 
quenched case as a function of the momentum q 2 . 
The line represents a fit of a power law to the data 
at q 2 < q 2 resulting in the exponent K w 0.2. 
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Figure 2: The (quenched) gluon propagator Das a 
function of q 2 . 



Figure 3: The gluon dressing function Z for the 
quenched and unquenched case as a function of q 2 
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Figure 4: The upper panel shows the zero three- 
momentum propagator C(f) of the gluon fields in 
quenched QCD at j3 =6.0 for a 48 4 lattice as a func- 
tion of time. In the lower panel the same data are 
shown, however, as logC(r) for C(t) > 0. 




Figure 5: The same as in the left figure, however, 
for the gluon propagator in full QCD on a 24 3 x 48 
lattice for one particular mass ma. 



the ghost and gluon propagators to our results using symmetric lattices (e.g. 48 4 ) we found large 
systematic effects at low momentum due to the asymmetry involved, in particular, for the lowest 
on-axis momenta along the elongated 'time' direction. 2 It is not excluded that sensible results at 
lower momenta for the gluon and ghost propagators can be extracted from asymmetric lattices by 
careful extrapolations. Attempts in that directions are made e.g. in [10]. 



3. Confinement criterion I: Violation of reflection positivity 

In view of the general objectives of this project it is also important to check whether our data 
for the propagators satisfy necessary criteria for confinement. For example, the gluon propagator 
has to violate reflection positivity, because otherwise gluons could be interpreted in terms of stable 
particle states. 3 To check whether reflection positivity is violated by the gluon propagator we have 
calculated the space-time (lattice) correlator at zero spatial momentum as function of time 

C « : =4= iW^expj^) (3.1) 

using our data for the gluon propagator D{Q,k^). The result for the quenched case (/? = 6.0, 48 4 
lattice) is shown in Fig. |4[ Obviously, reflection positivity is violated in a finite range of t. The 
same holds in the unquenched case as can be seen in Fig. ||, even though the lattice size in this case 
is too small such that the data do not bend over towards zero at larger t as in the quenched case 
(Fig. Q). Note that we have seen this happening also in our quenched data at jS = 6.0 using a 32 4 
lattice (not shown). 

2 For more details about this but also other systematic effects have a look at the Ph.D. thesis |Q| of one of us. 
3 Note that the ghost propagator violates reflection positivity trivially and hence ghosts are explicitly unphysical. 



4 



Lattice study of QCD Green's functions 



4. Confinement criteria II: The Kugo-Ojima confinement parameter 



Another criterion for confinement for QCD in covariant gauges was given by Kugo and Ojima 
long time ago [[Ij]. According to their scenario colored asymptotic states, if any, cannot be detected 
in the physical subspace due to the quartet mechanism. This is realized if a function u{p 2 ), defined 
as @,0] 



J d 4 xe^ ^D";c e (x)gf b " l A d v (y)c c (y) 



8^ 



PuPv 



8 ab u(p 2 



(4.1) 



has the zero-momentum limit 

-2\ 



U := lim u(p 



1 . 



(4.2) 



In Landau gauge this limit is connected to an in- 
frared diverging ghost dressing function J through 
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(4.3) 

Here v(p ) is an unknown function of p 2 (see [Q] 
for a definition). We have made an attempt to con- 
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firm the realization of the limit in Eq. ( |4.2[ ), not 
only by giving numerical evidence for a diverging 
ghost dressing function (see above), but also by 
estimating u(p 2 ) itself at different momenta p 2 in 
our lattice simulations 4 . Our estimates of u(p 2 ), 
renormalized at /I = 4 GeV, are shown in Fig. |6] as a function of p 2 . There we also show data of 
the ghost dressing function (renormalized at the same /i) in the form of an asymptote u defined as 

1 



i 

p 2 [GeV 2 ] 

Figure 6: Data of u(p 2 ,/j. 2 ) at j3 = 5.8 and 6.0 are 
shown using squares. Additionally, data of u are 
shown at the same j3 values (circles). All data re- 
fer to the same quenched configurations on a 32 4 
lattice and are renormalized at jj. = 4 GeV as de- 
scribed in [0| . Lines are drawn to guide the eye. 



u{P ,^ 



J(p 2 ^ 2 



1 



(4.4) 



Obviously, u(p ) and u(p ) are different at finite momentum, but according to Eq. ( |4.3| ) u(p ) has 



to reach asymptotically u(p 2 ) in the zero momentum limit. From our data shown in Fig. g we can 
confirm that the difference \u(p 2 ) — u{p 2 ) \ diminishes with decreasing momentum, even though an 
extrapolation of the given data to vanishing momentum is difficult to perform. For this, the explicit 
momentum dependence of u has to be known. Since u{p 2 ) seems to continuously approach u{p 2 ) 
with decreasing momentum and, by definition, u{p 2 ) is minus one at vanishing momentum for a 
diverging ghost dressing function our studied momentum range of u at least does not exclude the 



expected zero momentum limit (see Eq. (4.2)). 



5. The running coupling constant 

Given the dressing functions of gluon and ghost propagators a running coupling constant 
aM 2 ) = a^ 2 )Z(q 2 ,n 2 )J 2 (q 2 ,n 2 ) (5.1) 

4 For details on our estimation and renormalization of u(p 2 ) on the lattice we refer to Ref. [ph. 
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Figure 7: The running coupling oc s (q 2 ) as a func- 
tion of the momentum q 2 . 



based on the ghost-gluon- vertex can be calculated. 
This definition relies on the assumption that in 
Landau gauge this vertex stays bare also beyond 
perturbation theory. We have given numerical ev- 
idence in quenched and unquenched SU (3) lattice 
gauge theory confirming that the vertex renormal- 
ization constant Z\ « 1 in a (asymmetric) MOM 
scheme where the gluon momentum equals zero 
(see ^ 12]). Similar results indicating this, di- 
rectly and indirectly, were presented in lattice stud- 
ies of quenched SU (2) gauge theory [O], but also 
within the DSE approach [14]. Thus the definition 
of a s (q 2 ) according to Eq. (5J) is valid at least in 
that scheme. 

In Fig. we show our current data of OC s in the quenched case. One clearly sees it increasing 
towards lower momenta as long as q 2 > 0.3 GeV 2 . After passing a maximum at q 2 m 0.3 GeV 2 it 
decreases. We have found the same behavior on our sets of dynamical gauge configurations (see 
e.g. Fig. 4 in Ref. [12]). Therefore, on the basis of the present data we cannot confirm a s (q 2 ) to 
approach a non-trivial infrared fix-point, even less monotonously from below as expected from 
DSE studies. The reason for this is still unclear, however, in the light of the DSE results on a torus 
this might be a finite volume effect which we are unable to resolve at the present stage. 

Apart from the low-momentum region it is also interesting to look at a s {q 2 ) at larger momenta 
and to fit the corresponding 1-loop and 2-loop expressions to the data. This is shown in Fig. || for 
the quenched (l.h.s.) and unquenched case (r.h.s.). Obviously, the data follow the one and two- 
loop expression in the marked interval. However, large discretization errors are visible. Therefore, 
simulations at smaller lattice spacings are worthwhile in order to get more reliable values of A. 




q 2 [GeV 2 ] q 2 [GeV 2 ] 



Figure 8: The running coupling a s at large momenta for the quenched (left) and unquenched case (right). 
Since discretization errors are quite large we used only data from the 48 4 (24 3 x 48) lattice at /3 = 6.0 
(j3 = 5.29, ma = 0.0138) to fit the 1-loop and 2-loop expressions of a s to the data on the left (right) hand 
side of this figure. 
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6. Conclusions 

We have studied different aspects of lattice Landau gauge theory using MC simulations of 
quenched and full QCD. Thereby, we have mainly focused on the momentum dependence of the 
gluon and ghost propagators. We have demonstrated here for the first time that our lattice results for 
both propagators are consistent with different criteria for confinement, even though we could not 
confirm the infrared behavior as anticipated from DSE studies. Whether the gluon propagator at 
zero momentum vanishes in the limit of infinite volume or not has to be clarified in future studies. 

We could demonstrate that (dynamical) clover-improved Wilson fermions affect the gluon 
but only negligibly the ghost propagator. Importantly, the influence of fermions on the infrared 
behavior of both propagators seems to be small. 

On the basis of our present data for a s we cannot confirm it to approach a non-trivial infrared 
fix-point monotonously from below. In any case, our data at larger momenta agree with the 1-loop 
and 2-loop expressions of a s , even though simulations at larger /3 are worthwhile to perform to 
suppress discretization effects and to obtain more reliable results for A. 

All simulations were performed on the IBM pSeries 690 at HLRN and on the MVS-15000BM at the 
Joint Supercomputer Center (JSCC) in Moscow. This work was supported by the DFG under the contract 
FOR 465 (Forschergruppe Lattice Hadron Phenomenology), by the DFG-funded graduate school GK 271 
and with joint grants DFG 436 RUS 1 13/866/0 and RFBR 06-02-04014. We thank the QCDSF collaboration 
for providing us their unquenched configurations which we could access in the framework of the 13 Hadron- 
Physics initiative (EU contract RII3-CT-2004-506078). We are grateful to Hinnerk Stiiben for contributing 
parts of the program code. A. St. acknowledges discussions with Lorenz von Smekal. 
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